Determining the Dust Extinction of Gamma-ray Burst Host 
Galaxies: A Direct Method Based on Optical and X-ray 

Photometry 

Y. U\ Aigen Li^-^, and D.M. Wei^-^ 
ABSTRACT 

The dust extinction of gamma-ray bursts (GRBs) host galaxies, containing 
important clues to the nature of GRB progenitors and crucial for dereddening, 
is still poorly known. Here we propose a straightforward method to determine 
the extinction of GRB host galaxies by comparing the observed optical spectra 
to the intrinsic ones extrapolated from the X-ray spectra. The rationale for this 
method is from the standard fireball model: if the optical flux decay index equals 
to that of the X-ray flux, then there is no break frequency between the optical 
and X-ray bands, therefore we can derive the intrinsic optical flux from the X-ray 
spectra. We apply this method to three GRBs of which the optical and X-ray 
fluxes have the same decay indices and another one with inferred cooling break 
frequency, and obtain the rest-frame extinction curves of their host galaxies. The 
derived extinction curves are gray and do not resemble any extinction curves 
of local galaxies (e.g. the Milk Way, the Small/Large Magellanic Clouds, or 
nearby starburst galaxies). The amount of extinction is rather large (with visual 
extinction Ay ~ 1.6-3.4 mag). We model the derived extinction curves in terms 
of the silicate-graphite interstellar grain model. As expected from the "gray" 
nature of the derived extinction curve, the dust size distribution is skewed to 
large grains. We determine, for the flrst time, the local dust-to-gas ratios of 
GRB host galaxies using the model-derived dust parameters and the hydrogen 
column densities determined from X-ray absorptions. 
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1. Introduction 

It is widely acknowledged that the long-duration gamma-ray burst (GRB) is associated 
with the collapse of a massive star (Woosley 1993). Observational evidence supporting this 
coUapsar model includes the underlying supernova components in the afterglow of many 
GRBs (Zeh et al. 2004) and the observed location of GRBs in star-forming galaxies and 
active star-forming regions within their host galaxies (Paczynski 1998, Pruchter et al. 2006). 
In this scenario, GRBs are born and explode inside dense, dusty environments. The huge 
7-ray energy emission of GRBs is almost unaffected by absorptions, allowing them to be 
detected up to rather high redshifts (e.g. see Tagliaferri et al. 2005). Therefore, to study 
the dust and gas properties in the surrounding vicinity of GRBs is of great significance in 
understanding the interstellar medium (ISM) of star-forming galaxies throughout cosmic 
history. In addition, an accurate apprehension of the dust and gas immediate surrounding 
GRBs can also help (1) to reveal the nature of so-called "dark bursts" (i.e., whether the non- 
detection of some optical afterglow is due to dust extinction or the afterglow is intrinsically 
dark; see Lazzati et al. 2002 and references therein), (2) to detect the dust evolution with 
cosmic time, and (3) to correct for the extinction of optical emission in GRB afterglow 
analysis. 

The dust extinction of GRB host galaxies is traditionally modeled using either the Milky 
Way (MW), the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), or other 
presumed extinction curves (e.g. see Stratta et al. 2004; Kann et al. 2006; Schady et al. 2007; 
Starling et al. 2007; Tagliaferri et al. 2007). Recently, Chen et al. (2006) made the first effort 
to determine the extinction curves for GRB host galaxies without a priori assumption of 
the extinction law. The derived extinction curves differ from any known extinction laws of 
the Milky Way and external galaxies, challenging the traditional method commonly used in 
determining the extinction curves of GRB host galaxies. 

In this work we propose a novel, straightforward method to determine the extinction of 
GRB host galaxies by comparing the observed optical spectra to the intrinsic ones extrap- 
olated from the X-ray spectra. That such an analysis is possible follows from the standard 
fireball model. Based on the multi-wavelength afterglow photometry (including both the 
X-ray and optical data), we obtain the extinction curves of four selected bursts. We then 
model the size distribution and composition of the dust with the silicate-graphite interstellar 
grain model and obtain the dust-to-gas ratios in the local environment of GRBs. 
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2. Method 

The standard fireball model (Sari et al. 1998), which has been successful in explaining 
the overall properties of GRB afterglow (Mcszaaros & Rees 1997), predicts that the afterglow 
emission is produced by synchrotron radiation of electrons accelerated by the forward shock. 
In this model, with typical parameters, the optical to X-ray spectra can be described by a 
broken power law with indices (5 — {p — 1) /2 ior i/ < Uc or (3 — p/2 for v > i/c, where 
i/c is the cooling frequency and p is the electron energy distribution index. In most cases, 
the cooling break position is hard to determine. If the decay indices a of X-ray and optical 
bands are different, the cooling frequency lies between them, making it difficult to calculate 
the intrinsic optical flux from X-ray data. However, if the decay indices a of X-ray and 
optical bands are the same, then the optical and X-ray should lie on the same spectral 
segment, rendering it possible to calculate the intrinsic flux density in any optical band from 
Fx = Fx {\/Xx)^~'^, where /3 is the X-ray afterglow spectral index that we get from fitting 
the X-ray spectrum, and Fx is the X-ray flux density. After corrected for the Galactic 
extinction using the reddening maps of Schlegel et al. (1998), the observed spectral energy 
distribution (SED) of GRB at redshift z can be described as -Fa(i+z) = exp (— ^4^/1.086). 
Therefore, the extinction of the GRB host galaxy can be given by 

A, = 1.086 In ^^(VAx)^-^ _ 

With Ay interpolated, we can obtain the extinction curves (normalized to V band) of the 
GRB host galaxies. 

We then fit the derived extinction curve Ax/Ay with the standard silicate-graphite in- 
terstellar dust model which has successfully reproduced the extinction and IR emission of 
the MW galaxy, SMC and LMC (Weingartner & Draine 2001; Li & Draine 2001, 2002). 
The grain size distribution for both silicate and graphite is modeled with: dn = N (a) da <x 
a^""' exp(— a/ ac) da, where a is the grain radius (assumed spherical), ranging from Omin = 0.005 /xm 
to amax = 2.5/im, ac is the cut-off size. Note that it is assumed that both silicate dust and 
graphitic dust have the same size distribution. Let /gra be the number fraction of graphitic 
dust, the mass fraction of graphitic dust is /g^^ = /graPgra/ [/graPgra + (1 - /gra) Psii], where 
Psii ~ 3.5gcm~^ is the mass density of silicate material and pgra ~ 2.24 gcm~^ is that of 
graphite. 

With the fitted dust parameters, we can estimate the dust-to-gas ratio in each of the 
GRB host galaxies: 

rridust _ Mgra + Msii 
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where Nu is the hydrogen column density in the host galaxy; /iH is the atomic weight of H; 
the factor "1.4" accounts for helium; Mg^a and M^n are the column mass density of graphite 
and silicate material, respectively: 

Mgra = Nd -IT a^N (a) Pgra/gra dtt] (3) 

^ Ctmin 
/* ^max ^ 

= Na -n a'^N (a) psii (1 - /gra) (4) 



N (a) is the normalized dust size distribution; The dust column density A^^^ can be derived 
from 



Qrnax 

2 



Ax = 1.086 / N{a)TTa [/gra Qoxt,gra(a, A) + (1 - /gra) Qext,sil(a, A)] dtt Nd, (5) 

where (5ext,gra(o., A) and (5ext,sii(o, A) is the extinction efficiency of dust of radius a at wave- 
length A for graphite and silicate material, respectively. 



3. Data 

We select four GRBs that have both optical and X-ray observations. Photometric data 
are taken from literature (see Tabled 1,2). The optical to X-ray spectra are extracted when 
the afterglow light-curve are in a steady power-law state (e.g. see Panaitescu & Kumar 
2001, Fan Sz Piran 2006 for detailed analysis) to avoid complex phases (i.e. X-ray flares or 
re-brightening when the optical and X-ray emission are probably due to different components 
[Zhang et al. 2006]; see Fig. 1). For GRB 020405, GRB 030227 and GRB 060729, we adopt 
the spectra obtained when the cooling frequency Uc falls below the optical band, indicating 
an intrinsic single power law spectrum through optical and X-ray bands as discussed above. 
The decay indices a are all taken from literature except for GRB 060729 (around 0.35 day 
during plateau phase; for which a is not available in literature) we derive by fitting the 
afterglow light curve between 0.2-0.6 day. For GRB 061126, the decay indices of X-ray and 
optical bands are different, indicating a break frequency lying between them. At ~ 30 ks, 
the R band afterglow shows a break (see Fig. Id), which can be interpreted as the spectral 
break passing through the R bandjl] allowing us to calculate the intrinsic optical flux from 



^We note, however, that this burst, hke many other Swift bursts, does not obey the closure relation in the 
standard afterglow model (Perley et al. 2007), which adds uncertainties to our analysis. But the uncertainty 
of break frequency does not appear to affect the shape of the derived extinction curve - as can be seen in 
Figure 2, the extinction curve for the other three bursts remains gray even if we ignore the GRB 061126 
data. 
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4. Results 

We present in Table 3 the derived A\ of the GRB host galaxies at every observed optical 
band and in Table 4 the y-band extinction versus hydrogen column density {Ay/Nn) and 
the dust-to-gas ratio. The errors of the X-ray spectrum can bring about uncertainties on the 
extrapolated optical fluxes and thus on the derived ^4;^. We estimate the errors of Ax from P 
(the X-ray spectral index) through Eq.l. Larger errors of the X-ray spectrum result in larger 
uncertainties in Ax (e.g. see Table 3 and Fig. 2, GRB 020405). Most noticeably, the derived 
extinction curves of the four bursts are rather "gray" (see Fig. 2). Since all these extinction 
curves have very similar slope, we put all the extinction data of the four bursts together to fit 
them to the silicate-graphite grain model. The best fit parameters are 77 ~ 2.60, ac ~ 2.0 ^m, 

and ^ 0, with x' = EaiiAEaiiGRBs[(AMv)^od-(^AMi/)obs]Vc^' « 0.12 (obtained 
by summing up all wavebands and all GRBs, where a is the uncertainty for a given GRB at 
a given band). A prominent feature is the considerably small r] (the canonical value of 77 is 
3.5), indicating a grain size distribution skewed towards substantially large grains. The main 
reason for /gra ~ is the absence of the 2175 A extinction bump in the derived extinction 
curves which is generally attributed to small graphitic grains or PAHs. 

5. Discussion 

Previous works concerning dust extinction of GRB host galaxies mostly focused on 
fitting the observed photometry with the intrinsic power-low spectrum reddened by certain 
"standard" extinction curves inferred from the Milky Way or nearby galaxies (e.g. see 
Starling et al. 2007). However, lacking a priori knowledge of the dust properties in high 
redshift galaxies harboring GRBs. we have no reason to assume that they are the same as in 
local universe (e.g. see Stratta et al. 2007). Chen et al. (2006), for the first time, derived the 
extinction curves of GRBs without a priori assumption of the extinction law, but they only 
used the optical data. In this work, with carefully selected afterglow data covering X-ray to 
optical/near-infrared bands, we obtain the extinction curve of four GRB host galaxies more 
directly and precisely, based only on the standard fireball model. 

The "collapsar" model predicts GRBs to occur in active star-forming regions similar 
to Galactic molecular clouds (Jakobsson et al. 2006) which are heavily enshrouded by dust 
(Trentham et al. 2002; Tanvir et al. 2004). A recent dust scattering model proposed to 
account for the shallow-decay phase in Swift X-ray afterglow also requires large quantities 
of dust surrounding the GRBs (Shao & Dai 2007). Observations supporting the existence of 
large amount of dust include the emission and absorption features in some X-ray afterglows 
(Antonelli et al. 2000; Piro et al. 1999), large column densities of heavy elements revealed 
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by optical spectroscopy studies (Savaglio & Fall 2004; Savaglio 2006), and the non-detection 
at optical wavelengths for more than half the well- localized GRBs (Jakobsson et al. 2004). 
In contradiction with these evidence, traditional SED fitting often finds small extinction, 
primarily because the best fit model in most cases is the SMC-type extinction which, with 
a steep rise into the far ultraviolet (UV), often requires a small Ay to fit the spectrum 
(e.g. see Kann et al. 2006; Schady et al. 2007; Tagliaferri et al. 2007). Our work, showing 
considerably large Ay compared to that fitted with traditional method, is more consistent 
with theoretical prediction and observations. In addition, Rol et al. (2007) found that for 
GRB 051022 a lower hmit of Ay ~ 4.4 mag was needed, which imphes that at least in some 
GRBs the extinction Ay is rather large. 

The extinction curve derived in our work is flat, almost independent of wavelength, 
and is even "grayer" than the gray type of extinction curve obtained by Chen et al. (2006) , 
similar to the Calzetti et al. (1994) law suitable for local starburst galaxies. This result is in 
good agreement with other works fitting the SEDs of these bursts (e.g. Stratta et al. 2005; 
A. Li et al. 2007, in preparation). In particular, Pcrlcy et al. (2007) found that for GRB 
061126 the extinction curve is gray. Gray extinction has also been observed in Galactic dense 
clouds (Cardelli et al. 1988) and in the circumnuclear region of some AGNs (see Li 2007 for a 
review) . Gray extinction is produced by a dust distribution biased towards large grains (see 
§4), which may form from (1) grain coagulation naturally expected in the dense circumstance 
near GRBs (Maiolino et al. 2001a,b), (2) the biased evaporation of smaller grains due to the 
intense X-ray and UV radiation up to ~20 parsecs from the GRB (Waxman & Draine 2000, 
Pruchter et al. 2001, Savagho et al. 2003), and (3) preferential destruction of small grains by 
high energy ions in fast shocks (Jones 2004). Pcrna et al. (2003) computed the extinction 
curve that is obtained if standard Galactic dust is exposed to a GRB lasting more than a 
few tens of seconds (three of the four bursts in our sample meet this requirement, sec Tgg 
in Table 1) and found that the extinction curve can be very fiat, chiming with our result. 
We favor the grain growth hypothesis since the preferential destruction of small grains only 
occurs in the immediate GRB environment (~10-20 pc from the burst). 

It has long been proposed that GRB afterglow radiation, as well as the prompt emission, 
can destroy dust grains and cause Ay to decrease with time (e.g. see Vreeswijk 1999). We test 
this effect for GRB 060729, which is exceptionally bright in X-rays as well as at UV/optical 
wavelengths showing an unusually long unanimous plateau phase (~ Iday). We derive 
Ay 1.70 ± 0.20 mag Bit t — 0.35 day (in the plateau phase) and Ay ^ 1.59 ± 0.20 mag at 
t — 4.6 days (in normal decay phase) respectively (see Tables 1, 4), showing no significant 
dust destruction during this time. Detailed studies of dust destruction by GRBs will be 
presented in a forthcoming paper (Z. Jin et al. in preparation). 
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In accordance with previous works, we find that the average Av/A^h is smaller than that 
m the Milky Way@ which is usually ascribed to a lower dust-to-gas ratio in GRB vicinities 
(e.g. see Watson et al. 2006). However, there is no obvious reason why the amount of dust is 
low in the dense environment surrounding GRBs. We note that the dust extinction is very 
sensitive to the dust size distribution, for larger grains the extinction (on a per unit mass 
basis) is low, but the amount of dust may be still high (e.g. see Li 2007). In fact, based on 
the model fit dust parameters, the dust-to-gas ratios for most bursts are larger than that in 
the Milky Way. On the other hand, grain growth through coagulation in dense molecular 
clouds enshrouding GRBs is expected and this would result in a dust size distribution biased 
in favour of large grains, a flat extinction curve, and a reduced Ay /Ny{. 
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grants 10621303 and 10673034, and the National Basic Research Program of China (973 
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Fig. 1. — Light curves of the four selected bursts. The first three panels (a, b and c) are bursts with the 
same decay index in both X-ray and optical bands. The last panel (d) is for GRB 061126, for which the 
X-ray and optical bands have different decay indices and cooling frequency between them. Vertical lines 
denote the time when the adopted GRB afterglow spectra were obtained. 



- 12 - 




2 4 6 8 10 

MX (|im) 



Fig. 2. — Rest-frame extinction curves (normalized to the V band) for 4 GRB host galaxies and the best-fit 
extinction curve from the silicate-graphite model (thick black line). Also shown are the Milky Way, SMC 
and LMC extinction curves. 
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Table 1: Observational Properties of the four GRBs 



GRB 


z 


(i 


Fx (A^Jy) 


t (days) 


Tgo (s) 




ax 


References 


020405 


0.691 


1.0 ±0.2 


0.23 


1.98 


60 


1.54 ±0.06 


1.97 ± 1.10 


1,2 


030227 


4 


0.94 ±0.05 


0.125 


0.87 


18 


0.95 ±0.16 


0.97 ±0.07 


3,4 


060729 


0.54 


1.06 ±0.01 


0.2 


4.6 


115 


1.27 ±0.10 


1.29 ±0.03 


5 


061126 


1.1588 


0.5 ±0.07 


40.5 


0.023 


191 


0.75 ±0.06 


1.31 ±0.01 


6 


060729*1 


0.54 


1.06 ±0.01 


3.46 


0.35 


115 


0.26 ±0.07 


0.35 ±0.15 


5 



Note. — z is the redshift of the burst; /3 is the intrinsic optical/UV to near-lR spectral index derived from 
the standard afterglow model; Fx is the X-ray flux density at 1 keV; t is measured from the burst trigger 
time; Tgo is the duration of the burst; and ax, respectively the temporal decay index at optical and X-ray 
bands, are all taken from literature except GRB 060729, for which we obtain ao and ax by fitting the the 
afterglow light curve between 0.2-0.6 day. 

References. — (1) Berger et al. 2003; (2) Stratta et al. 2005; (3) Castro-Tirado et al. 2003; (4) Mereghetti 
et al. 2003; (5) Dirk Grupe et al. 2007; (6) Perley et al. 2007. 
Data are taken during the plateau phase. See §5 for discussion 



Table 2: Optical/UV to near-IR flux density 



GRB 












(AiJy) 












UVM2 


UVWl 


V 


B 


V 


R 


/ J 


H 


K 


020405 








5.93 


7.18 


9.5 


10.9 


15.7 28.6 


34.9 


42.5 


030227 










2.4 




2.7 




10.2 


15 


060729 


7.56 


10.18 


11.05 


21.98 


23 


25.79 










061126 




65.12 


54.66 


97.41 


120.18 


143.18 


173.6 


217.97 309.38 


409.42 


523.16 


060729* 


124.2 


170.1 


195.3 


313.1 


359.4 


455.1 











Note. — The flux densities F^ measured in the observer frame are taken at the time of the vertical lines 
reported in Fig. 1. All the data have been corrected for Galactic extinction 



Tabic 3: Derived extinction in the observed wavelength A. 



GRB 
























UVW2 


UVM2 


UVWl 


u 


B 


V 


R 




.7 


H 


K 


020405 








2.66 ± 1.12 


2.65 ± 1.16 


2.60 ± 1.21 


2.70 ± 1.26 


2.58 ± 1.31 


2.29 ± 1.39 


2.37 ± 1.45 


2.47 ± 


030227 










2.88 ± 0.26 




3.16 ± 0.28 






2.58 ± 0.33 


2.54 ± 


060729 


1.91 ± 0.05 


1.71 ± 0.05 


1.82 ± 0.05 


1.41 ± 0.05 


1.62 ± 0.06 


1.71 ± 0.06 












()(ill2(i 




4.f)') ± (i.o:! 


4.:i8 ± f).():i 


:',-<).'! ± 0.04 


'■•,.>>\ ± 0.04 


:!.74 ± 0.04 


:i.fi:! ± 0.04 


:i.40 ± 0.04 


:i.:',4 ± o.o.'i 


:',.!<) ± 0.0.1 


:',.()(! ± 


()()IJ72!)"- 


2.U2 ± O.lJ.'i 


1.7!) ± O.U.'i 


i.hi = u.o.rj 


i.(i7 ± 0.0."") 


i.78 ± 0.(J.rj 


i.74 I (J.O.'i 
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Table 4: V^-band extinction and dust-to-gas ratio in the rest frame of each burst. 



GRB 


Av(mag) 


AThCIO^^ cm-^) 


Av /Nuii-O^'-''-' magcm^) 


'^dust/"'gas(10 


020405 
030227 
060729 
061126 


2.60 ± 1.17 

2.57 ± 0.32 
1.59 ± 0.20 
3.35 ± 0.43 


0.8 ± 0.2 

0.076 ± 0.003 
1.1 ± 0.3 


3.2 ± 1.6 
O.asto'ii 
20.9 ± 2.8 

3.0 ± 0.9 


0.99 ± 0.62 

IQ + O-O'' 
"■-^■^-0.03 

6.5 ± 0.9 

0.94 ± 0.28 


060729* 


1.70 ± 0.20 


0.076 ± 0.003 


22.4 ± 2.8 


6.9 ± 0.9 



Note. — Av is the rest-frame V-band extinction of GRB host galaxies; A^h is the rest-frame equivalent 
column densities of hydrogen measured from X-ray absorption assuming a solar metal abundance at the 
same time when the multi-band spectra were taken (see §2) 



